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Abstract. The analytical transfer matrix technique is applied to the Schrodinger 



equation of symmetric quartic-well potential problem in the form V(x) = ^kx 2 + Ax 4 



This gives quantization condition from which we can calculate the ground-state energy 
eigenvalues numerically. We also compare the results with those obtained from 
numerical shooting method, perturbation theory, and WKB method. 
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Figure 1. The single-stepped potential 



1. Introduction 

Quantum Mechanical bound-state problems have long been of interest to physicists. 
There exist several means to study them, e.g. WKB approximation [8], time- 
independent perturbation theory [8], numerical shooting method [H], finite element 
method [12, 13J. Here we analyze one-dimensional problem using the analytical transfer 
matrix method (ATMM) devised by pQ. This method originated from planar optical 
wave-guides [2] , and tunneling [3] . The main principle is to divide the domain into many 
tiny segments. Each segment possesses a constant potential. The concept of a transfer 
matrix then arises when we connect the wave functions at the boundary of two different 
potential levels. 

2. Single-stepped potential 

For a background , we consider the stepped potential with a particular energy eigenvalue 
as shown in Fig.l. It is well-known that the wave functions are plane waves. The original 
idea stems from the fact that the wave functions and their first derivatives can generally 
be written as linear combination of the form 
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where a, /3, 7, r/, are constants to be determined. We may substitute plane wave solution 
into 
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Here we regard Mj 



a /3 
7 ^ 



as the matrix that transfers wave function and its first 



derivative from the position x = Xj-i to x = Xj. Solving Eq.(2) and Eq.(3) altogether, 
we obtain 



a = cos(Kj_ici), (3 



sin(« 3 -_id), 



7 = sin(/c 3 -_id), 77 = cos(/«j_id). 
Substituting Eq.(4) into Eq.(l), we have 

cos(/«j_id) sin(«j_id) 

sin(«j_i(i) cos(/«j_id) 



V'i-i(^i-i) 
V^-ife-i) 
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where 



v /2m(E - V£ 



j = 1,2,... 



On the other hand, in case the energy is lower than the stepped potential we would 
rather get 



V>j-l(Xj-l) 

^•-iOj-i) 



cosh(«j_id) ^-j- sinh(«j_ 1 d) 
sinh(/Cj_id) cosh(/Cj_id) 
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3. The multi-stepped potential well 



In order to solve the Schroedinger equation in each of three regions in Fig. 2(a), the 
boundary conditions due to the continuity of the wave function and its first derivative 
are applied at the boundaries of regions I, II, and III. In case the energy is less than 
both V\ and V3, the solution is given by [H] 

ka = mix — sin" 1 [P 2 ] — sin -1 [Pi], (7) 

where 
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Figure 2. (a) Asymmetric finite square- well potential in case of no phase contribution, 
(b) The stepped potential well with phase contribution. 



Consider the stepped potential well shown in Fig. 2 (b). We pay attention to the case 
of E > Vi. The wave function has the from of Aoexp (Pqx) in the regime x < xq and 
/^exp (—P 3 x) in the regime x > d\ + d 2 . Between them, the wave function oscillates 
with different wavelengths whenever the particle moves from a constant potential to 
another step. We finally obtain the quantization rule as follows pQ: 



K\d\ + fi^G^ + &( s ) — n7T + tan 





+ tan 1 
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(8) 



where 



$(s) = $ 2 — tan" 

The second and third terms on the right-hand side of Eq.(8)are half-phase losses 
at the potential barriers V and V 3 , respectively. We observe that by setting V\ = V 2 , 
we obtain $(s) = 0, this phase contribution actually results from the interference of the 
scattered sub- waves between the potentials V\ and V 2 . 



K2 
Kl 



tan($ 5 



nir + tan 
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4. An arbitrary potential-well function 

Since a continuous potential well may be viewed as a stack of thin films each of which 
possesses a constant potential, the above method can therefore be extended to study an 
arbitrary one- dimensional potential well of the form shown in Fig. 3 below: 

We have also truncated the profile at x = xq and x = x s . The transaction certainly 
affects the values of energy levels as compared to the situation in idealized system. The 
effects will clearly be negligible if the potential at the transaction points is very much 
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Figure 4. A graph of an arbitrary potential well is equivalent to an assembling of tiny 
bars. 



larger than energies of relevant levels. Given that x c and Xd are turning points, we 
divide the region (xo, x c ),(x c , xj) and (xd, x s ) into I, m and n equal parts each of width 
d. According to the classically allowed or forbidden regions, the multi-stepped potential 
corresponding to the ith, jth and fcth section layers of the form are shown in Fig. H] 
The transfer matrices can be written as 

cosh(aid) — i sinh(ajd) 
— ati sinh(aj<i) cosh(aj(i) 



Mi 



12 J- 



(9) 



M,- 



cos (Kjd) — — sin( Kj d) 
Kj sin ( Kj d) cos ( Kj d) 



,j = I + 1,1 + 2,..., I + m; 



(10) 
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and 



where 



a. 



cosh(afc(i) — sinh(o!fc(i) 
— a>k smh(akd) cosh(afcci) 



y/2m[y(x t ) - E] 

h 



k = i +771+I, l+m+n; (11) 



y/2m[E - V(xj)] 



h 



(12) 



^2m[V{x k ) - E] 



h 



Applying the boundary conditions at (x = x ) and (x = x s ) yields 

I _ l+m _ l+m+n 



ip(x ) 

i)'(xo) 



n M - n«, n m 



i=i 



j=l+l 



k=l+m+l 



ip(x s 
i[)'(x s 



(13) 



where the prime denotes differentiation with respect to x. As known, the wave functions 
in the region x < Xq is A e p °^ x ~ x °^ and wave function in the region x > x s is A s e~ Ps<yX ~ Xs ' > 
where 



P 
P, 



^2m[V 


~E] 


h 


y/2m[V s 
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X = Xq, 



(14) 



Aq, A s are the amplitude coefficients to be determined. From Eq. (12) the solution for 
j = I + m by using similar algebra manipulation as developed in the reference is given 

by m 

'Pl+m+1 



Ki +m d = ni +m -K + tan 



-i 



- K 2+r 
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Summing all indices j, we have 

l+m 

Kjd+$(s) = A/V+tan" 1 

j=l+l 
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where 



l+m-l 

$( s ) = J^+i-tan" 1 

3=1+1 



K j+ i 
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(15) 



A" = 0,1,... (16) 



(17) 



and 



1), = tan 1 (— ), Pj = /v, tan 



tan 



3+1 



K/jd 



$(s) is the phase contribution devoted by the scattered sub-wave. p and p +m+ i 
are equivalent exponential decaying coefficients corresponding to the regions x < x c 
and x > Xd, respectively, Pi is momentum that transfers energy at a stack of thin films 
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interval i - 
interval k 



1,2, ... , / and is momentum that transfers energy at a stack of thin films 
: l + rn+l,l + m + 2,...,l + m + n. We now investigate the half-phase losses 



and tan 



l + m+l 



in Eq. (16). It is 
at X[ + i and Ki +r 



at the turning point, i.e., the two terms tan 1 

obviously clear that as l,m, and n — > oo, d — > 0, we have = at xi + i and Kz +m = 
at xi +m . Pi and P; +m+ i are positive and finite for bound states so that the half-phase 
losses at the turning point have the value |. We thus obtain a quantization condition 
as the width of the section layers d tends to zero, i.e. 



K(x)dx + $(s) = (J\f+ 1)7T, TV 



0,1,2, 



5. Numerical Results 



We now pay attention to the particular problem of calculating ground-state energy of 
the potential-well 

V(x) = ~mcu 2 x 2 + Ax 4 (19) 

where A is a positive constant. After introducing new variable and parameter as 
(3 = , £ = ax, a = (^fO 5 , Schrodinger equation is transformed to 

^l+[e-e-(3em) = o. (20) 

Also, the potential in term of new variable and parameter is given by 

v(o = e+Pt 4 - (21) 

Due to this, the principal phase in the quantization condition Eq. (18) reduces, for the 
sake of computation, to 

I d n{x)dx = f d Vz-e- pe dt, (22) 

J X c J £ c 

where £ c = (^) 5 x c , = {^)^ x d- The ATMM gives a quantization rule similar to 
that obtained from WKB method, the only difference is that for ATMM the phase 
contribution can be calculated with slight complication whereas WKB result may be 
equivalent to the ATMM case by setting the phase contribution equal to ^. We have 
calculated the ground-state energy eigenvalues under this potential as shown in Table 
1. 



6. Conclusion 



The NSM is generally regarded as one of the most efficient methods that give accurate 
results because it integrates the Schrodinger equation directly, though in numerical 
sense. From tablel, we see that the ATMM we have adopted from Cao et al pQ 
gives outstandingly better results for ground state than those obtained from l si order 
perturbation theory and the typical WKB method for every value of A (or (3) under 
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Table 1. Comparison of ground-state energy obtained from ATMM., standard WKB., 
I s * order perturbation, and numerical shooting method. 



3 = ATMM WKB l st P.T. 5 NSM 5 

(ft*) /2) (ftv/2) (ftj/2) (ftj/2) 



0.1 


1 


.06640625 


1 


.03515625 


1.07500000 


1. 


.06451896 


0.2 


1 


.12131250 


1 


.06656250 


1.15000000 


1 


.11740450 


0.3 


1 


.16906250 


1 


.09515625 


1.22500000 


1 


.16304919 


0.4 


1 


.21234375 


1 


.12171875 


1.30000000 


1 


.20371079 


0.5 


1 


.25156250 


1 


.14687550 


1.37500000 


1 


.24065919 


0.6 


1 


.28781250 


1 


,16934375 


1.45000000 


1 


.27469829 


0.7 


1 


.32171875 


1 


,19106875 


1.52500000 


1 


.30637709 


0.8 


1 


.35390625 


1 


.21203125 


1.60000000 


1 


.33609065 


0.9 


1 


.38359375 


1 


.23187550 


1.67500000 


1 


.36413529 


1.0 


1 


.41183595 


1 


.25156250 


1.75000000 


1 


.39073967 



this quartic single well potential. In our viewpoint, the complication of the ATMM 
computation lies at the evaluation of phase contribution, particularly when the potential 
is not of very simple form so that we cannot obtain analytical expression from the 
quantization rule but numerical one instead, in which we employ iteration method to 
work them out. In fact, the ATMM is claimed by Cao et al pQ to give exact formalism 
of quantization rule without any approximation. One slight disagreement is emphasized 
here in that the truncation of the bound potential at xq,x s in fact indicates a kind of 
approximation. Anyhow, this estimation seems to be a reasonable one that does not 
considerably affect the energy calculation especially when we choose x s far away enough 
from the RHS turning point. 
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